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We investigate the equilibrium properties of a quasi-two-dimensional degenerate boson-fermion 
mixture (DBFM) with a bosonic vortex-antivortex superposed state (VAVSS) using a quantum- 
hydrodynamic model. We show that, depending on the choice of parameters, the DBFM with 
a VAVSS can exhibit rich phase structures. For repulsive boson-fermion (BF) interaction, the 
Bose-Einstein condensate (BEC) may constitute a petal-shaped "core" inside the honeycomb-like 
fermionic component, or a ring-shaped joint "sheU" around the onion-like fermionic cloud, or muhiple 
segregated "islands" embedded in the disc-shaped Fermi gas. For attractive BF interaction just 
below the threshold for collapse, an almost complete mixing between the bosonic and fermionic 
components is formed, where the fermionic component tends to mimic a bosonic VAVSS. The 
influence of an anharmonic trap on the density distributions of the DBFM with a bosonic VAVSS 
is discussed. In addition, a stabihty region for different cases of DBFM (without vortex, with a 
bosonic vortex, and with a bosonic VAVSS) with specific parameters is given. 

PACS numbers: 03.75.Lm, 03.75.Ss 



I. INTRODUCTION 

With the experimental realizations of atomic Bose- 
Einstein condensate (BEC), molecular BEC and 
fermionic condensate, considerable attention has been 
given to degenerate boson-fermion mixtures (DBFMs) 
where particles obeying different quantum statistics are 
intermingled ^^. The DBFMs, such as '^Li-^Li [IJl, 
23Na-6Li g, 87Rb-40K [^, and "4Yb-i73Ybllll, 
have been successfully observed by different experimen- 
tal groups via interspecies sympathetic cooling. In ad- 
dition, such mixtures can also be realized from an im- 
balanced two-component fermi gas where all minority 
fermions pair up with majority fermions and form a BEC 

MM- 

The coUisional interactions between bosons and 
fermions strongly affect the properties of DBFMs, even 
though the ultracold mixed gases are very dilute J14| - t31| . 
A noticeable example is the structure and stability of a 
trapped DBFM which has been investigated extensively 
by several authors. In the case of repulsive boson-fermion 
(BF) interaction, the system is expected to undergo a 
mixing-demixing transition with the increase of the re- 
pulsion, while in the case of attractive BF interaction 
the mixture collapses above a critical strength [y, [3, [iJ- 
[l9l . |33 - |35{ . In particular, the effect of a bosonic vortex 
on the stability and the mixing- demixing transition of a 
DBFM is recently addressed [23[ ISg]. It is shown that 
when the BEC sustains a vortex state the DBFM be- 
comes more stable. 

In this paper, we consider the equilibrium proper- 
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ties of a quasi-two-dimensional (quasi-2D) DBFM in the 
presence of a bosonic vortex-antivortex superposed state 
(VAVSS). Recently, the VAVSS in BECs has attracted 
extensive research interest as it exhibits peculiar petal- 
like structure and rich dynamics |37h41| . Furthermore, 
the VAVSS may have fundamental as well as practical 
applications in quantum information and inertial sens- 
ing. Experimentally, the creation of VAVSS in BECs has 
been reported by Anderson et al [42| and Wright et al 
l43l| . Throughout the present work, we use a quantum- 
hydrodynamic model [3^ to investigate the structure of 
the exact 2D spatial density distributions and the sta- 
bility of the system against collapse. This mean-field- 
hydrodynamic model is quite successful in the study of 
collapse dynamics [Ig] and dark 44 1 and bright soli- 
tons [H, Hi] in a DBFM. The theoretical predictions 
on fermionic collapse in a DBFM agree well with the 
experimental results (y, |9| , and those on bright solitons 
are in agreement with a relevant microscopic investiga- 
tion \4(y\. We show that, for interspecies repulsion, the 
bosonic component may form a petal-shaped "core" in- 
side the honeycomb-like fermionic one, or a ring-shaped 
joint "shell" around the onion-like fermionic cloud, or 
even segregated "islands" embedded in the disc-shaped 
Fermi gas. For interspecies attraction below the critical 
value for collapse, there is an enhanced mixing (even an 
almost complete mixing) between bosons and fermions, 
where the fermionic cloud tends to simulate a bosonic 
VAVSS. Moreover, we find that the DBFM with a bosonic 
VAVSS and the one with a bosonic pure vortex sustain 
an intermittent stability diagram. 

The paper is organized as follows. In section II we in- 
troduce the mean-field-hydrodynamic model for a quasi- 
2D DBFM, which is comprised of two coupled quantum- 
hydrodynamic equations. In section III, we present the 



numerical results and some discussion on the structure of 
the exact 2D density profiles of a DBFM with a bosonic 
VAVSS. Furthermore, a stability region for different cases 
of DBFM with specific parameters is displayed. The con- 
clusion is outlined in the last section. 



II. MEAN-FIELD-HYDRODYNAMIC MODEL 

FOR A DEGENERATE BOSON-FERMION 

MIXTURE 

We consider a dilute degenerate mixture composed of 
Nb condensed bosons and Np spin-polarized fermions at 
zero temperature, with respective particle mass ttib and 
mp. In the case of strong axial (z-direction) confinement, 
the DBFM reduces to a quasi-2D system, where the axial 
trap frequency is denoted by Wz. Here we use the time- 
dependent mean-field-hydrodynamic model developed in 
Ref. [3a] to describe this BF system. The Lagrangian 
density of the DBFM reads 



C = Cb+Cf+C 
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Here Cb is the ordinary bosonic Lagrangian density. 
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where ^B(a^,2/,i) is the hydrodynamic field of the Bose 
gas, i.e., the macroscopic BEC wave function, ub = 
{"^B I is the 2D bosonic density with normalization Nb = 
//I^bI dxdy, and Vb is the external trapping poten- 
tial for the bosons, qbb = 2\/2Trh^aBB/{azB''nB) is 
the 2D interatomic interaction strength with ubb being 
the 3D boson-boson (BB) s-wave scattering length and 
QzB — ^/^/i'^B^z) the axial harmonic length of bosons. 
The bosonic Lagrangian density Cb can describe well all 
the d yna mical properties of the dilute quasi-2D BEC (see 
Ref. (36| and references therein). For simplicity, here we 
have assumed the BEC to be strictly 2D. 

The Lagrangian density of the Fermi gas Cp is ex- 
pressed by [3a] 



Cp = 
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where \I'i?(x, j/,i) is the hydrodynamic field of the de- 
generate Fermi gas (DFG) (i.e., an average fermionic 
wave function), np = ]^f1 is the 2D fermionic den- 
sity with normalization Np — Jf \'i'p\ dxdy, and Vp is 
the external trapping potential for the fermions. £,p{np) 
denotes the zero-temperature bulk energy density of an 



ideal quasi-2D Fermi gas under axial harmonic confine- 
ment, and it reads 



hujzTr{npazp)'^, 
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(4) 
For ^ n^- < Xji^ixa^^p), the Fermi gas is strictly 2D, 

where a^p — ^Jh/(rapujz) is the axial harmonic length 
of fermions. For np ^ l/(27ra^p), the Fermi gas has 
the 2D-3D crossover, where several single-particle modes 
of the harmonic oscillator along the z axis are occupied. 
The Fermi gas becomes 3D for np ^ l/(27ra^^). Due 
to the Pauli exclusion principle, the interaction between 
identical fermions in spin polarized state is highly sup- 
pressed and has been neglected in the Lagrangian density 
Cp and will be neglected throughout this paper. The sec- 
ond term — h^ \^x,y'^F\ /Qnip in Eq.Q represents the 
Thomas-Fermi- Weizsacker kinetic energy which includes 
the surface kinetic energy due to spatial variation [l9| . 
Note that the fermionic Lagrangian density Cp can only 
be safely used to describe the static and collective prop- 
erties of the quasi-2D Fermi gas [3g. As the strict 2D 
Fermi gas requiring a very small number of particles is 
not the case of real experiments, we take into account 
the 2D-3D crossover in the fermionic component. 

It may not be entirely proper to define an average 
fermionic wave function '^ p for a DFG. A proper treat- 
ment of the DFG should be performed using a fully anti- 
symmetrized many-body Slater determinant wave func- 
tion. However, the probability density np of a DFG 
calculated in this way should lead to reasonable results 
IJj and has led to proper probability distribution for a 
DFG [l6|, uM as well as results for collapse of a DFG 



[la . Il9l l32j ] in agreement with experiment. This method 
has also been used successfully to predict fermionic bright 
and dark solitons in a DBFM [2^ i44] • The virtue of the 
mean-field-hydrodynamic model for a DFG over a mi- 
croscopic description is its simplicity and good predic- 
tive power. Similar mean-field treatment for degenerate 
fermions can also be found in Refs. (47l.l48|. 

Finally, the Lagrangian density Cbp for BF interaction 
is given by 



Cbf — 



'QBFnBnp, 



(5) 



where gBP — 2'KaBpfi^/{'mBP\/'2TrazBO'zF) is the cou- 
pling constant of BF interaction in terms of the 3D BF 
s-wave scattering length ubf and the BF reduced mass 
niBP = mBmp/{mB + nip). 

In terms of Eqs. (IT])-([5]), the Euler-Lagrange equations 
of motion become 
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Here ^xp = d^p /dnp is the bulk chemical potential of the 
non-interacting Fermi gas in the 2D-3D crossover, which 
is expressed by 
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The resultant quantum-hydrodynamic equations have a 
nonpolynomial nonlinearity for the fermionic component 
as shown in Eqs. dH), and dH). 

For simplicity, we assume that the 2D radial trapping 
potentials for the bosonic and fermionic components are 
the same and take the forms 



1 
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VBix,y) = Vp{x,y) = -ttibuj 

(9) 

where uq = y^h/{mBUJ±) is the radial harmonic length of 
bosons with uj±_ the radial trap frequency as seen by the 
bosons. The quartic term in Eq.(j9l) denotes the anhar- 
monicity of the traps, which is possibly caused by exper- 
imental uncertainties or artificial quartic confinement. A 
is a dimensionless parameter that determines the relative 
strength of the quartic component for bosons. Our choice 
of the trapping potentials corresponds to a reduction of 
w_L by a factor y/mB/rnp as in the study by Modugno et 
al [34I and a reduction of A by a factor ^JrriBjmp, where 
the two assumptions give a simpler analytical form of the 
final dynamic equations without any consequence to our 
qualitative study. Recent studies show that the quartic 
distortion influences the rotational properties [49| and 
the collective excitations [50| of BECs. An interesting 
question is whether the additional quartic distortion in- 
fluences the structure and stability of a DBFM. 

Eqs. (6) and (7) can be renormalized to dimension- 
less form. By introducing the notations xq = x/a^^ 
2/0 = y/ao, V = tuj±, p = niB/mp, c = QbbNb/ {huj±al), 
f = gBFNp/[hujj_al), k = gBENslifu^^al), ipB = 
ao'^B/y/NB and V-F = o-o^ p/^/Np, we obtained the 
rescaled Euler-Lagrange equations for the mixture 
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where rijo = IV'jl {j — B,F) are the rescaled bosonic 
and fermionic densities. Here the rescaled bulk chemical 
potential of the Fermi gas with 7 = ujz/(^± reads 
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For the case of a DBFM sustaining a bosonic vortex, a 
mixing-demixing transition and a collapse in the mixture 



have been predicted by Adhikari et al \3u\ . In view of the 
recent studies of VAVSS in BECs ,37.-,4a] it is of inter- 
est to see how the structure and the stability of a DBFM 
modify when the Bose gas is in a VAVSS. The appearance 
of a quantized VAVSS as well as a pure vortex state is 
the genuine confirmation of phase coherence and super- 
fiuidity of the boson component in the DBFM. As the 
VAVSS can stably exist in nonrotated BECs [33, [111, we 
consider the bosonic VAVSS of the dimensionless form 

EI 

i^B = Aer(--o+y'o)/^-'e'S^a{xo + tyoY + ^e'H^o - iyoY], 

(13) 
where A is the normalization constant, a denotes the 
width of BEC, (5 is a constant phase which may be taken 
to be zero without loss of generality, and I is the winding 
number describing the quantum circulation. The real 
parameters a, /? characterize the proportion of the vortex 
and antivortex with a'^ + /3'^ — 1. The relative phase (f> 
between the vortex and the antivortex just causes offset 
of the density distribution by an angle (j)/2l. 



III. RESULTS AND DISCUSSION 

Here we numerically solve the coupled quantum- 
hydrodynamic equations (jlOp and (jlip with bulk chem- 
ical potential p,Qp given by Eq. (J12p and a trial bosonic 
VAVSS with definite ratio of a^ and /3^. The trial so- 
lution of the fermionic component is chosen as a Gaus- 
sian wave function. In view of the petal-like structure 
of the VAVSS in BECs, we perform exact 2D numerical 
calculations which require an enormous computation ef- 
fort. In our numerical simulation, we consider the ^^Rb- 
^'^K mixture with tub being the mass of ^^Rb atom and 
nip the mass of ^°K atom, and take u}z/2Tr = 100 Hz, 
7 — 10. In order to investigate the structure and stabil- 
ity of the DBFM with a bosonic VAVSS, we seek the sta- 
tionary solution of Eqs. (ITU)) and pT|) in the presence of 
a bosonic VAVSS (fT5)) . There are several typical numer- 
ical approaches for obtaining the stationary solution of a 
nonlinear Schodinger equation, for instance, the Newton 
relaxation method [51[and the imaginary time propa- 
gation (ITP) method '38', '53|. In this paper, we obtain 
the ground state of the DBFM with a bosonic VAVSS by 
using the ITP method based on the split-step Fourier al- 
gorithm [53| . Our preliminary study shows that the petal 
structure of a VAVSS with a^ 7^ /3^ only emerges under 
the condition of extremely weak interatomic interaction 
or very small particle number, which is usually not met in 
experiments. For a VAVSS with a^ — /S"^ — 1/2, however, 
the petal structure always appears for arbitrary value of 
the interatomic interaction or the particle number. So we 
shall mainly consider the case of a^ — /3'^ — 1/2. The rel- 
ative phase is taken to be (/) = throughout this paper. 
The BB and BF interactions can be arbitrarily tuned 
by the Feshbach resonance technique through varying a 
background magnetic field [2|. 



Firstly, we consider the solutions of system in the ab- 
sence of interspecies interaction. In figure 1 we show 
the 2D bosonic density distribution nBoixo,yo) (top), the 
2D fcrmionic one npQ(xo,yo) (middle), and the column 
density distributions njo{xo,yo — 0) (j — B,F) of the 
bosonic and fermionic components (bottom) in a DBFM 
with a bosonic VAVSS, where Nb = 1000, Np = 100, 
Z = 1, a^ = /3^ = 1/2, QBE = 40 nm, asF = 0, and 
the anharmonic parameters are A = (left) and A = 0.1 
(right) , respectively. The bosonic component in the mix- 
ture exhibits a petal-like structure (see figures 1(a) and 
(d)), which is similar to the case of a pure BEC. At 
the same time the fermionic component displays a Gaus- 
sian distribution due to the absence of BF interaction, as 
shown in figures 1(b) and (e). Comparing with the case 
of a harmonic potential (A = 0, see figures 1(a) and (b)), 
the density distributions nBoixQ^yo) and npQ(xQ,yQ) in 
the presence of an anharmonic term (A = 0.1, see figures 
1(d) and (e)), i.e., in the presence of a quartic distor- 
tion, become more tightly confined. This can also be 
seen in the column density distributions tijq {xq ,yo ~ 0) 
{j — B, F) along the xq axis (figures 1(c) and (f)). Thus 
the additional quartic distortion may favor the sympa- 
thetic cooling between the bosons and the fermions, and 
makes the mixture more stable. Our numerical simula- 
tion shows that when A 7^ 0, 0.1 the density distributions 
"■50(2^07^0) a-iid npo{xo,yo) £^re similar to those in the 
case of A = 0, i.e., the basic structure of the DBFM is 
not affected by the quartic distortion. The conclusion is 
the same in the presence of BF interaction. So in the 
following we illustrate the equilibrium properties of the 
DBFM with a bosonic VAVSS with A = 0. 

When the BF interaction is nonzero, the DBFM with a 
bosonic VAVSS can display rich phase structures. Figure 
2 gives the density distribution npo{xo,yo) of fermions 
(left) and the column density distributions njo{xo,yo = 
0) {j ~ B. F) of bosons and fermions (right) for various 
BF interactions, where A = 0, A^s = 1000, Np = 100, 
I = \t a^ — P^ ~ 1/2, aBB = 40 nm. As shown in figures 
2(a) and (d), for repulsive aBF = 90 nm the system is 
in a partial mixed state, where the bosonic component 
with a petal-like structure lies inside the fermionic one 
while the fermionic cloud shows a honeycomb-like struc- 
ture due to the BF repulsion. For a sufficiently large re- 
pulsive aBF = 190 nm, the fermionic cloud is completely 
expelled outside the BEC (see figures 2(b) and (e)), which 
is usually referred as demixing, i.e., the system is in 
a complete phase-separated state. The phenomenon of 
mixing-demixing transition has also been found in previ- 
ous theoretical studies of a DBFM JJ,ia,Ml or a DBFM 
with a bosonic vortex [3^ . For the DBFM with a bosonic 
VAVSS, however, with the increase of BF repulsion the 
DEC can not be completely expelled from the central 
region of the harmonic trap due to the bosonic VAVSS. 
In the case of attractive BF interaction [aBP < 0) we 
find that the fermionic component is pulled inside the 
bosonic one and the maximum mixing between the two 
components is then achieved for a critical value of ubf- 
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FIG. 1: (color online) Bosonic density distribution UBoi^o, J/o) 
(top), fermionic density distribution nFoi^o,yo) (middle), 
and column density distributions njo{xo,yo — 0) {j = B,F) 
(bottom) for a DBFM with a bosonic VAVSS in the absence 
of BF interaction. The anharmonic parameters are A = 
(left) and A = 0.1 (right), respectively. Here Nb = 1000, 
Nf = 100, I = 1, a'^ = P^ = 1/2, and ubb = 40 nm. 



With the further increase of the strength of the attrac- 
tive BF interaction the system collapses. In figures 2(c) 
and (f), where ubf — ^90 nm, just below the thresh- 
old for collapse, we can see an almost complete mix- 
ing between the bosonic and fermionic clouds. In this 
mixed state, the fermionic density distribution developes 
a petal-like structure reminiscent of a VAVSS as in the 
bosonic component. However, near the center of the trap 
the fermionic density npo{xo = 0,yo) tends to a nonzero 
constant value and does not have the VAVSS behavior. 
For a VAVSS with I ^ 1 anda^ = ^'^ = 1/2, the probabil- 
ity density in the gap regions between two petals should 
be zero. Therefore the fermionic component is not really 
in a VAVSS but tends to mimic a bosonic VAVSS due to 
mixing. 

Depending on the choice of parameters the separated 
phase may exhibit different configurations. Shown in fig- 
ure 3 are the bosonic density distribution nBo{xo,yo) 
(top), the fermionic one npo{xo,yo) (middle), and the 
column density distributions 'njo{xo,yo = 0) (j = B,F) 
(bottom) in a DBFM with a bosonic VAVSS for the case 
of equal particle numbers Nb = Np = 1000. Here A = 0, 
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FIG. 2: (color online) Fermionic density distribu- 
tion nFo{xo,yo) (left) and column density distributions 
njo{xo,yo = 0) (j = B, F) (right) of a DBFM in the presence 
of a bosonic VAVSS, where A = 0, A^s = 1000, Nf = 100, 
; = 1, a^ = ^^ = 1/2, and ubb = 40 nm. The BF scattering 
lengths are ubf = 90 nm (top), ubf = 190 nm (middle), and 
GBF = —90 nm (bottom), respectively. 



Z = 1, a^ = /32 = 1/2, GBB = 2000 nm, and the BF 
s-wave scattering lengths are asp ~ 500 nm (left) and 
aBF — 2000 nm (right), respectively. It is easy to con- 
firm that the diluteness condition of the gases is satisfied 
for all the parameters mentioned in this paper. From 
figures 3(a)-(c), we can see for aBF = 500 nm and given 
parameters the system is in a perfect mixed state, where 
the fermionic density reaches the maximum in the gap 
region between two bosonic petals because of the BF re- 
pulsion. For sufficiently strong BF repulsion asp — 2000 
nm, the bosons are completely expelled from the center 
of the trap as shown in figure 3(d), forming a ring-shaped 
joint "shell" around the fermions due to the presence of 
the VAVSS. At the same time the fermionic component 
develops an onion-like core, which can be seen in figure 
3(e). The column density distributions njQ{xo,yo — 0) 
(j = B,F) are illustrated in figure 3(f). Figures 3(d)-(f) 
definitely indicate a new separated phase of the DBFM 
which is different from that in figures 2(b) and (e), where 
though the fermions are completely expelled outside the 
bosons there is remarkable fermionic density along the 
axis of a^o = in the trap center. 

In figure 4, we display the bosonic and fermionic den- 
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FIG. 3: (color online) Bosonic density distribution nBo{xo, j/o) 
(top), fermionic density distribution nFo(2:o,J/o) (middle), 
and column density distributions njo{xQ,yo = 0) {j = B,F) 
(bottom) for a DBFM with a bosonic VAVSS in the case of 
Nb = Nf = 1000. The BF s-wave scattering lengths are 
aBF ~ 500 nm (left) and ubf = 2000 nm (right), respec- 
tively. The other parameters are I — 1, a'^ — jS^ — 1/2, and 
aBB = 2000 nm. 



sity distributions njo(xo,2/o) {j = B,F) as well as the 
column density distributions njo{xo,yo — 0) (j = B,F) 
for the case of Nb = 100 and Np — 1000, where the 
interspecies s-wave scattering lengths are aBP — 800 nm 
(left) and Ubf — 8000 nm (right), respectively. The 
other parameters are X — 0, I — 1, a^ = [3^ = 1/2, 
and qbb = 800 nm. For repulsive qbf = 800 nm the 
DBFM begins to evolve into a coexistence phase of com- 
ponent mixing and separation, where the bosonic compo- 
nent forms two distant and segregated "islands" embed- 
ded in the disc-shaped fermionic gas (see figures 4(a)-(c)). 
When gbf = 8000 nm the two bosonic "islands" become 
more compact and closer, and the fermionic component 
is completely separated from the bosonic one (see figures 
4(d)-(f)). As shown in figure 4(f), the fermionic com- 
ponent constitutes both a "shell" around and a "core" 
inside the BEC along the xo-direction, i.e., the Fermi 
gas distributes on both sides of the Bose gas along the 
XQ-direction, which is similar to that found in Ref. jl4|. 
However, the layer distribution of the present separated 



phase is only along the xo-direction in the nearby re- 
gion oi yo — and is absent along the j/o-direction, i.e., 
the present separated phase does not show 2D or 3D 
layer configuration. Furthermore, the core structure or 
the core plus shell structure of the fermionic cloud in a 
DBFM without vortex can be achieved by varying the BF 
interaction continuously [14], whereas the different sepa- 
rated phases in the DBFM with a bosonic VAVSS corre- 
spond respectively the different parameter conditions and 
can not be achieved by merely change the BF coupling 
continuously. The above analysis indicates that the inlay 
configuration of the two bosonic "islands" embedded in 
the disc-shaped Fermi gas is a novel separated phase. We 
expect that this novel separated phase can be observed 
and tested in the future experiments. For large asF the 
fermions could possibly form p-wave pairing, which is 
beyond the scope of the mean- field treatment. The two 
compact bosonic segregated "islands" in figure 4(d) per- 
haps can be understood as a consequence of buoyancy in 
the Fermi sea. 

The BF interaction has a very strong influence not only 
on the structure of the DBFM with a bosonic VAVSS 
but also on the stability of the system. For a repulsive 
BF interaction (positive ubf), a stable structure can al- 
ways be achieved for a fixed Nb, Np, X, I, a, 13, ip, and 
BB scattering length ubb as shown in figures 2-4. How- 
ever, for a sufficiently strong attractive BF interaction 
(negative qbf), the system can undergo a simultane- 
ous collapse of the bosonic and fermionic density distri- 
butions. Physically, the critical value of BF scattering 
length is determined by the balance between the repul- 
sion of bosons and fermions and the mutual attractive 
BF interaction. Likewise, when the boson number Nb or 
fermion number Np are sufficiently large, the attractive 
BF interaction can not be stabilized by the repulsions 
in BB and BF subsystems. Thus the mixture lowers its 
energy via increasing the boson and fermion densities, 
and finally the bosonic component or the fermionic one 
or both the bosonic and fermionic components collapse 
simultaneously due to instability. In our investigation, 
the instability signature is found by monitoring the fail- 
ure of the numerical iterative process, i.e., the invalida- 
tion of the coupled equations ^TU\\ and (ITT]) describing 
the DBFM with conserved quantities Nb and Np. In 
fact, an indefinite growth of the maximum of the bosonic 
and fermionic densities always indicates the appearance 
of an instability onset, and therefore the collapse of the 
mixture is triggered and unavoidable. This approach is 
accurate and reliable |23j . 

The stability region in the plane spanned by \og{Np) 
and \og{NB) is shown in figure 5. As an illustration, we 
assume the BB and BF scattering lengths of the ^^Rb- 
^"K mixture are qbb = 5 nm and asp = — 15 nm which 
are very close to the typical experimental values gbb = 
5.25 nm and qbf — —13.8 nm reported by Ferrari et al 
|54| . The four different curves (square, circle, triangle, 
and inverted triangle) in the figure 5 mark the stability 
limit for four different cases of DBFM: (i) DBFM without 




FIG. 4: (color online) Bosonic density distribution nBoi^o, yo) 
(top), fermionic density distribution nFo{xo,yo) (middle), 
and column density distributions njo{xo,yo ~ 0) {j = B,F) 
(bottom) for a DBFM with a bosonic VAVSS in the case of 
Nb = 100 and Nf = 1000- The BF s-wave scattering lengths 
are asF = 800 nm (left) and asF ~ 8000 nm (right), respec- 
tively. The other parameters are I — 1, a'^ — (3^ — 1/2, and 
ciBB ~ 800 nm. 



vortex, (ii) DBFM with a bosonic VAVSS oil = 1 and 
^2^/32^ 1/2, (in) DBFM with a bosonic VAVSS of 
I :^ 1, a^ ^ 3/4 and /S^ = 1/4, and (iv) DBFM with a 
bosonic vortex of Z = 1, respectively. The region below 
each curve denotes the stability region of the mixture, 
and that above each curve corresponds to the instability 
(collapse) one of the system. The stability curves in the 
range of \og{Np) ^ 2 are four almost overlapped and 
parallel horizontal lines (data not shown here) which are 
similar to those in the range of log(iVi7') ^ 3.9. 

From figure 5 we can see for the same parameters the 
DBFM with a bosonic vortex is more stable than the one 
without vortex due to the repulsive centrifugal kinetic en- 
ergy, in agreement with previous study [36] . In addition, 
for the case of Np s^ int(103-«) « 6310 the DBFM with 
a bosonic vortex is more stable than that with a bosonic 
VAVSS because the energy of a VAVSS (especially, a 
VAVSS with equal proportions |a| = |/3|) is higher than 
that of a pure vortex (or antivortex) for a given / [3q . 



However, for the case of 6310 < iV^ < int(103-9) w 7943 
we find that the DBFM with a bosonic VAVSS is more 
stable than that with a bosonic vortex. Notice that for 
attractive BF interaction the fermionic density distribu- 
tion tends to simulate the bosonic one (see figure 2(c)), 
thus the overlap region between the bosons and fermions 
in a DBFM with a bosonic VAVSS is smaller than that 
in a DBFM with a bosonic vortex due to the petal-like 
structure of VAVSS. This effect dominates the stability 
of the system when the number of fermions is sufficiently 
large such that the influence of the BF attraction on the 
collapse of a DBFM with a bosonic VAVSS is weaker than 
that of a DBFM with a bosonic vortex, which is a prob- 
able reason for the above intermittency phenomenon in 
the stability region. At the same time, there is also an in- 
termittency phenomenon between the stability curve for 
the DBFM with a bosonic VAVSS of a^ = 1/2 and that 
of a^ ^ 1/2 (e.g., a^ — 3/4). The analysis is similar to 
the above discussion. Note that the stability curve for 



that th e different separated phases in the DBFM w ith 



the DBFM with a bosonic VAVSS of ^ = 1, 



= 1/4 



and (3^ = 3/4 is the same with that of Z = 1, a^ = 3/4 
and /3^ = 1/4 because the two systems have the equal 
energy [38|, which is verified in our numerical simula- 
tion. Furthermore, the DBFM with a bosonic VAVSS 
(/ = 1) is more stable than the conventional one without 
vortex (/ — 0) in the wide range of Np < 7943 because 
for the bosonic VAVSS there is an additional repulsive 
centrifugal term fi^P /(2mB{x^ + y^)) in the stationary 
GP equation [38!|. When Np ^ 7943 all the four different 
cases of DBFM become instable for Nb above the order 
of unit. For the DBFMs with a bosonic vortex (VAVSS) 
of higher Z, our simulation shows that the stability curves 
are similar to the case of Z = 1 , but the systems of higher 
I are less stable than those oil = 1. 



IV. CONCLUSION 

We have investigated the structure and stability of a 
quasi-2D DBFM with a bosonic VAVSS, where the BB 
interaction is taken to be repulsive and the BF inter- 
action to be both repulsive and attractive. It is shown 
that the bosonic and fermionic density distributions for a 
DBFM with a bosonic VAVSS in a harmonic plus quartic 
potential are more tightly confined than those in a har- 
monic potential. The additional quartic distortion due 
to experimental uncertainties or artificial quartic confine- 
ment does not infiuence the stationary structure of the 
mixture, and may favor the sympathetic cooling between 
the interspecies and makes the system more stable. 

Depending on the choice of parameters, the DBFM 
with a bosonic VAVSS can show rich phase structures. 
For repulsive BF interaction, the bosonic component can 
form a petal-shaped "core" inside the honeycomb-like 
fermionic one, or a ring-shaped joint "shell" around the 
onion-like fermionic cloud, or multiple segregated "is- 
lands" embedded in the disc-shaped Fermi gas. Note 
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log(N^) 

FIG. 5: (color online) Stability region for the degenerate 
*^Rb-'"'K mixture with scattering lengths asB = 5 nm and 
aBF = — 15 nm. Here the curves of square, circle, triangle, 
and inverted triangle show the stability limit for the DBFM 
without vortex, the DBFM with a bosonic VAVSS oi I = I 
and c? = /3^ 



1/2, the DBFM with a bosonic VAVSS of 



; = 1, a^ = 3/4 and /3^ = 1/4, and the DBFM with a bosonic 
vortex of / = 1, respectively. For log(A'fl) above the limit or 
log(7VF) on the right of the limit the bosonic component or 
the fermionic one or both the two components collapse simul- 
taneously. 



a bosonic VAVSS are formed under different parameter 
conditions and can not be achieved by just changing the 
BF coupling continuously, which is different from the 
case of a conventional DBFM [ij]. In addition, there 
is a mixing-demixing transition which is controlled by 
the strength of the BF repulsion. For attractive BF 
interaction below the threshold for collapse, an almost 
complete mixing between the bosonic and fermionic com- 
ponents can be formed, where the fermionic component 
tends to simulate a bosonic VAVSS. Furthermore, we 
give the stability region for four different cases of DBFM 
with specific parameters, where the two stability curves 
of the DBFM with a bosonic VAVSS and the one with 
a bosonic vortex display an intermittent configuration. 
These new findings can be verified in the future experi- 
ments on DBFMs, and thus the present work provides a 
new way to test further the validity and prediction of the 
quantum-hydrodynamic model. 
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